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ABSTRACT On the basis of conclusions drawn from previous experimental and theoretical evidence, 
we extend the two-fluid solar-wind model by including the collisionless dissipation of 
hydromagnetic waves originating at the sun. We generate a series of solar wind models 
parameterized by the total energy flux of hydromagnetic waves at the base of the model. 
The resulting properties of propagation and dissipating of hydromagnetic waves on this 
model are presented. 


We interpret the strong positive correlation between 
observed solar wind speeds vg (subscript E refers to 
values at 1 AU) and proton temperatures T„g [Burlaga 
and Ogilvie, 1970; Hundhausen et al., 1970] as defining 
a continuum of average macroscopic states. In this case, 
a realistic model of the wind should be capable of 
predicting a continuous range of wind speeds and proton 
temperatures consistent with this Tpg-vg correlation. 
Since the wind is frequently observed to blow much 
faster than the “quiet day” wind, it convects much more 
energy at some times than at others. We believe that this 
variability in energy flux reflects the dissipation of 
varying amounts of nonthermal energy supplied to the 
system, a view consistent with the Tp£-vg correlation. 

On this basis, we have recently shown [Hartle and 
Barnes, 1970], in the context of the two-fluid model 
[Sturrock and Hartle, 1966; Hartle and Sturrock, 1968] 
that it is possible to choose, ad hoc, a class of 
hypothetical energy deposition functions that give 
models consistent with the Tpg-vg correlation. Ac- 
cordingly, we found that primary energy deposition 
occurs over an extended region up to about 25 R Q from 
the sun. This result is compatible with the nonthermal 
heating mechanism suggested by Barnes [1968, 1969] ; 


that is, the collisionless dissipation of fast-mode 
hydromagnetic waves originating at the sun. 

Here we extend our previous model by replacing the 
artificial heat source with that corresponding to 
hydromagnetic dissipation. We present some of the 
properties of propagation and dissipation of hydro- 
magnetic waves as found from a self-consistent 
treatment of the two-fluid equations with such heating. 

This treatment is simplified by assuming that the solar 
wind is composed of protons and electrons which 
undergo steady, radial, spherically symmetric expansion. 
The flow is then described in terms of the proton 
density n and flow speed v (equal to electron density 
and flow speed), the electron temperature T e , and the 
proton temperature Tp as functions of the heliocentric 
distance r. These profiles are determined by the 
two-fluid equations of continuity, momentum, and 
energy given in Hartle and Barnes [1970] by 


nvr 2 = J = constant 


nm p v %^~i [nk(T e +Tp)}~ 


GM e mpn 




-vkT e 


dn 

dr 



( 1 ) 

( 2 ) 


The authors are with the NASA Laboratory for Planetary 

Atmospheres, Goddard Space Flight Center, Greenbelt, Md., and = ” 2 l> E n ^'^e ~ 

the Space Science Division, Ames Research Center, Moffett 
Field, California. 


248 



to I'w 


.004 


dT„ 


nvk 


dr 


vkT. 


dn 


P dr 


1 d 
r 2 dr 



= J v E n k{T e ~T p )+ 


( 4 ) 


in terms of the proton and electron masses trip and m e , 
the solar mass M e , the gravitational constant G, and 
Boltzmann’s constant k. The electron and proton 
thermal conduction coefficients, K e and Kp, and the 
electron-proton energy exchange rate vg are those of 
Braginskii [1965], The fast-mode hydromagnetic heat 
sources are of the form 
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The local plasma-frame (complex) circular frequency, in 
terms of the rest-frame (of the sun) frequency co 0 and 
wave vector k, is given by co = or + i y = o> 0 -k • v + il 
and is related to k through the dispersion relation 
[Barnes, 1968]. The factor W is the energy density 
spectrum of the waves and G corresponds to the energy 
absorption by the protons (electrons) per unit wave 
energy in the time 1/co. The heat functions are based on 
the assumptions that the velocity distributions are 
maxwellian, that the wave amplitudes are small enough 
that linear wave theory is valid, and that (-yl oj. 

The energy deposition term p\ a complicated 

w 

function of the plasma parameters n, v, T e , and Tp, the 
local average magnetic field B (r), and the angle 6 
between the wave vector k and field B. We assume that 
the average properties of the magnetic field BfrJ are 
adequately described by the spiral model of Parker 
[1963]. In the models considered here, the electron 
energy deposition term G e =* 0 since the electron- 
electron collision frequency is much greater than the 
source frequency co 0 = 2X10" 2 sec" 1 (see below). The 
typical dependence of the proton energy deposition 
term Gp on 9 is shown in figure 1 for plasma parameters 
at the base of one of the models (see below; 
F 0 = 5.2X10 3 ergs cm" 2 sec" 1 for this case). Here we 
note that the energy deposition distribution peaks at 
9 ~20° and that primary proton heating occurs over the 
range 10° <0 <35°. 



Figure 1. Proton energy absorption rate per unit wave 
energy in time 1/co versus 6, the angle between the wave 
vector and magnetic field. 

Since the wavelengths in question are small relative to 
macroscopic scale lengths, we determine the energy 
density W of the waves in the approximation of 
“geometrical hydromagnetics.” In this case, the ray 
paths x(s) and wave vectors k( s) (s is path length along 
ray) are obtained from Hamilton’s equations, the 
dispersion relation, and the computed profiles n, v, Tp, 
T e , and B. Then for a given pair (x 0 , k 0 ) at the base and 
fixed c o Q 
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where 7^ = Im is the damping decrement and 3co/9k 
the group velocity, and the integration is along the ray. 

The solar-wind models considered here are determined 
by self-consistent numerical solutions of equations (1) 
through (7) and the ray equations using a method of 
iteration similar to that described in Hartle and Barnes 
[1970]. In addition to specifying the magnetic field 
strength B 0 and the hydromagnetic wave spectrum at 
the base of the model, the solutions are subject to the 
usual constraints of specifying the base density and 
temperatures n 0 , T eo , Tp 0 , requiring the velocity v to 
pass continuously through the subsonic-supersonic 
transition, and requiring T e and Tp to tend to zero as r 
goes to infinity. We treat the wave spectrum as a discrete 
number of rays (six in this work), weighted so that the 
base intensity is approximately isotropic in the outward 
directions (zero in the inward directions), and 
the frequency spectrum is assumed monotonic. 
The boundary values we choose at the base 
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r- 2R e are n Q = 1.5X10 6 cm" 3 , T eo = 1.3X10 6 °K, 
Tp 0 - 1.7X10 6 °K, and B 0 = 0.18, consistent with 
coronal observations [Newkirk, 1967]. Our selection 
Tpo T eo , implying selective proton heating below the 
base, is consistent with the inner corona heating model 
of D’Angelo [1968, 1969]. The source frequency is 
taken to be co o = 2X10 -2 sec -1 , near the peak of the 
observed photospherxc and chromospheric acoustic 
spectrum [Leighton et al., 1962; Tanenbaum et al., 
1969] . 

We have calculated a series of models that are 
parameterized by the remaining free parameter of the 
system; namely, the total energy flux F 0 of hydro- 
magnetic waves at the base. Three of the resulting ray 
trajectories are shown in figure 2 in terms of their 
azimuthal and radial coordinates 0 and rjR Q . 
These trajectories, corresponding to a flux 
F 0 = 5.2X 10 3 ergs cm - 2 sec -1 , are typical of the class of 
models we considered and, for purposes of illustration, 
have been selected so that they all become radial at 
0=0. The base emission angles d e = cos -1 (k • r/k), 
between the wave vectors and the radial direction, are 
0°, 15°, and 45°. We note that the ray with emission 
angle d e = 0 propagates radially while the remaining rays 
refract toward the radial direction, becoming essentially 
radial at 10 or 12 R & . In addition, we observe that the 
refraction becomes stronger as the emission angle 
increases. 


40 - 



Let us now consider the damping of these rays. The 
peak amplitude of the energy deposition term Gp 
increases strongly with ftp = 8nnkTp/B 2 ; for example, 
the damping rate y « exp(-l/j3p) when fip < 0.1. 
Nevertheless, since the shape of the energy deposition 
term Gp of figure 1 is essentially the same for the range 
of parameters of interest, we can get a qualitative 


indication of where the relative damping between rays 
takes place by considering the angles 8 each wave vector 
makes with the magnetic field. The 8 profiles 
corresponding to the rays of figure 2 are shown in figure 
3. In this case, the ray with emission angle 8 e = 45° 
should lose a significant portion of its energy at radii less 
than about 6 R & since it passes through the range 
10° < 8 < 35° where the damping is strongest (damping 
will be relatively weak for r < 2.5 R e ). The intermediate 
ray, 8 e = 15°, will not damp as strongly near the base, 
resulting in a more extended heating region. Finally, the 
radial ray corresponds to the most extended dissipation 
since it has relatively low values of 8 throughout. The 
rays dissipate at similar rates beyond about 12 R 0 where 
they are essentially radial. In this region 8 increases with 
the spiral angle, resulting in increasing deposition rates 
with r, here the remaining portion of the wave energy is 
lost. 
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Figure 3. Angle 8 between the wave vector and mag- 
netic field versus radial position r/R e . 

Upon summing the heat contribution from each ray 
we obtain the proton heat source Pp on this model. The 
resulting energy deposition distributions are shown in 
figure 4 corresponding to three models with 
F 0 = 2.2X10 3 , 5.2X10 3 , and 1.2X10 4 (ergscm -2 sec -1 ), 
which result in wind speeds vg = 290, 370, and 427 
(km sec -1 ) and proton temperatures Tpg = 2X 1 0 4 , 
6.2X10 4 , and 1.7X1 0 s (° K), respectively. These models 
are consistent with the T p g- Vg correlation as discussed 
earlier by Barnes and Hartle (p. 219). For clarity, we 
normalized the heat profiles by the density to give the 
heating rate per proton. In each model hydromagnetic 
heating is important over an extensive region of 
20-30 R 0 in radius. Beyond this region, Pp decays 
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Pp/n [ergs/sec] 



rapidly and other terms in the proton heat equation 
become dominant. 

When the energy flux F 0 is increased, we note in figure 
4 that there is an overall increase in the heat profile 
resulting in both increased wind speeds vg and proton 
temperatures Tpp at 1 AU. This is consistent with our 
previous model [Hartle and Bames, 1970], which 
demonstrated that heating in the region of subsonic flow 
(r < 5-7 R q ) primarily raises the wind speed at 1 AU and 
heating in the supersonic flow region primarily raises the 
proton temperature at 1 AU. Also consistent with a 
requirement of our previous model is the inward motion 
of the heat source peak with increasing flux F 0 . This 
variation in the peak is related to the fact that increased 
heating raises $ p - SmkTp/B 2 , which in turn raises the 
damping rate 7 <*exp(-l IPp), leading to a relatively 
larger rate of energy deposition near the base. 
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